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Abstract. We report the detection of faint emission in 
the high-excitation [Oiv] 25.90/im line in a number of 
starburst galaxies, from observations obtained with the 
Short Wavelength Spectrometer (SWS) on board ISO. 
Further observations of M 82 spatially resolve the [O iv] 
emitting region. Detection of this line in starbursts is sur- 
prising since it is not produced in measurable quantities 
in H II regions around hot main-sequence stars, the domi- 
nant energy source of starburst galaxies. We discuss vari- 
ous models for the formation of this line. [O iv] that is spa- 
tially resolved by ISO cannot originate in a weak AGN and 
must be due to very hot stars or ionizing shocks related 
to the starburst activity. For low-excitation starbursts like 
M 82, shocks are the most plausible source of [Oiv] emis- 
sion. 



Key words: Galaxies: ISM 
frared: galaxies 



Galaxies: starburst - In- 



1. Introduction 

Mid-infrared fine structure lines are powerful probes of 
dusty and obscured galactic nuclei, being able to pene- 
trate extinctions up to the equivalent of Ay ~ 50. Using 
the Short Wavelength Spectrometer (SWS) on board the 
Infrared Space Observatory (ISO), it is possible to detect 
faint lines and sources. The rich observed spectra can be 
used for a detailed modelling of the ionizing spectra of 
starbursts (e.g. Rigopoulou et al. 1996| , Kunze et al. 1996 ) 
and AGNs (Moorwood et al. 1996 ). Clear differences be- 
tween their spectra make these lines a valuable new tool 
for discriminating between AGN and starburst activity in 
visually obscured galaxies. AGN spectra include emission 
from highly ionized species and the so-called coronal lines, 
requiring photons up to ^300 eV for their creation. In con- 
trast, starburst spectra are dominated by lines of low ex- 



Send offprint requests to: D. Lutz (lutz@mpe.mpg.de) 

* Based on observations with ISO, an ESA project with in- 
struments funded by ESA Member States and with the par- 
ticipation of ISAS and NASA. The SWS is a joint project of 
SRON and MPE. 



citation species, because even hot, massive stars emit few 
ionizing photons beyond the Hen edge at 54 eV. 

Line ratios like [Oiv]25.9^m / [Nc n] 12.8/zm and 
[Nev]14.3/im / [Neil] 12.8/im have been used by Lutz 
et al. (1996a) and Genzel et al. (1998) to establish the 
dominant source of luminosity in ultraluminous infrared 
galaxies (ULIRGs). In some of the starburst templates 
studied, very faint [O iv] emission was found, about two 
orders of magnitude weaker than in typical AGNs. Faint 
[O iv] emission in starbursts is not relevant for establish- 
ing the power source of ULIRGs, but its origin poses an 
interesting problem because its creation ionization energy 
is slightly above the Hell edge. In this letter we examine 
possible mechanisms for its production. 

2. Observations and data reduction 

Observations of a variety of starburst galaxies in the 
[O iv] , [Ne il] and [Ne in] lines have been obtained with 
the ISO-SWS in 1996 and 1997, as part of a more com- 
prehensive guaranteed time program. In addition, we use 
data from a raster of SWS observations along the major 
axis of M 82 obtained on March 16, 1996 in an open time 
program. One of the target lines of this program was [Fe n] 
25.988/im, which is close enough in wavelength to extract 
the [O iv] line from the same scans. 

Integration times of our observations in SWS02 mode 
were typically 2 seconds per step, i.e. 200 seconds for 
the complete up-down scan covering the line. We used 
standard procedures from the SWS Interactive Analysis 
system for data reduction. For strong sources like bright 
starbursts, residual fringing is often obvious in the pro- 
cessed data. This was corrected for by fitting sine func- 
tions to line-free regions of the spectra. This procedure 
is sufficiently reliable for the brightest sources, since the 
fringes can be approximated by a single sine function over 
the small observed wavelength range near [O iv] , and the 
lines are narrow. Nevertheless, baseline uncertainty due 
to fringing is often the limiting factor in measuring the 
[O iv] line flux, and allows us to set only upper limits in 
some sources. [Ne n] and [Ne in] were always significantly 
stronger than residual fringes at these wavelengths. 
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Table 1. Observed line fluxes in starburst galaxies 



25.80 25.85 25.90 25.95 26.00 
Wavelength [/im] 

Fig. 1. [O iv] towards the southwest starburst lobe of M 82 
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Fig. 2. [Oiv] 25.9 /im emission relative to the lower exci- 
tation starburst emission for our sample galaxies, plotted 
against starburst excitation as given by the [Ne in] 15.5 /im 
/ [Neil] 12.8 /xm line ratio 



3. Results 

Fig. [j] presents the first detection of [O iv] emission in a 
starburst, in the spectrum of M 82. The line fluxes mea- 
sured for all our sources are summarized in Table [j], which 
also compares the [O iv] flux to the flux of lower excitation 
starburst lines. For this comparison, we use [Neil] 12.81/iin 
+ 0.44x [Nem] 15.55/im, the factor 0.44 chosen to give 
equal weight (by mass) to singly and doubly ionized neon. 
Contrary to normalizing to just [Ne n] or [Ne ill] , this mea- 
sure will be robust to changes in excitation of the star- 
burst proper. For two galaxies without measured neon 
lines (NGC 6764 and NGC 6052), we have estimated the 
neon fluxes from the measured [S in] 33.48/im flux using 
the average scaling for the other sample galaxies. 

Accurate linewidths are difficult to determine for the 
faint [Oiv] lines. All the detected lines are consistent in 
width with the starburst fine structure lines in the same 
source, that is unresolved or slightly resolved at the SWS 
resolving power of ^1000. 

Faint [O iv] emission is thus fairly universally detected 
in starburst galaxies, at the percent level compared to 
the starburst neon lines. There is a weak correlation of 
[O iv] strength with excitation of the starburst (Fig. |J) . 
High-excitation, low-metallicity starbursting dwarfs like 
II Zw 40 and NGC 5253 exhibit relatively stronger [O iv] 



Source 


F([0iv]) 
10" 20 Wcm" 2 


IOIVJ 
[NoII]+0.44x [NcIII] 


M 82 


8.00 


0.008 


II Zw 40 


0.55 


0.062 


NGC 253 


5.00 


0.012 


IC 342 


<1.00 


<0.011 


NGC 3256 


0.93 


0.012 


NGC 3690 A 


<1.20 


<0.041 


NGC 3690 B/C 
NGC 4038/39 1 


0.80 


0.024 


<0.90 


<0.087 


NGC 4945 


1.40 


0.017 


M 83 


0.80 


0.006 


NGC 5253 


0.65 


0.038 


NGC 6052 


<0.80 


<0.127 


NGC 6764 2 


<0.80 


<0.228 


NGC 7552 2 


<1.20 


<0.018 



1 Interaction region 

2 [Ne n] and [Ne ill] fluxes estimated from observed [S ill] (see 
text) 

than low-excitation starbursts, but the effect is not very 
pronounced compared to the large excitation difference as 
measured by [Ne in] /[Neil]. 

3.1. Spatially resolved [Ow] emission in M 82 

The SWS02 raster of M 82 was obtained along the major 
axis (PA 68°), using a spacing of 9". At the time of the 
observation, the position angle of the long axis of the SWS 
apertures was 60.9°, i.e. almost aligned to the major axis. 
We compared the measured [O iv] spatial profile (Fig. || 
top) to the SWS beam profile along the long dimension 
of the SWS aperture, for the relevant wavelength band 
(A. Salama, priv. communication). The [Oiv] emission is 
clearly resolved compared to the point source beam profile 
and must originate in a region similar to the size of the 
entire starburst region. Further support that [Oiv] is re- 
solved comes from the observation that its 'rotation curve' 
follows the rotation of M 82 (Fig. || bottom). 

4. Discussion 

At a level of just a percent of the strongest low excitation 
lines, a variety of possible excitation mechanisms for the 
[O iv] line must be considered. 

4.1. Weak AGNs 

Because of the great strength of [O iv] emission in Seyfert 
galaxies, quite faint and perhaps obscured AGNs embed- 
ded in a more luminous starburst would contribute suffi- 
cient [O iv] , with the additional constraint that their nar- 
row line width would have to be small. In fact, hard X-ray 
observations may indicate AGNs deeply hidden in some of 
our sources. The case is convincing for NGC 4945 (Iwa- 
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n=1000, Matthews and Ferland continuum, stop at T=100K 



20 
NE Offset [arcsec] SW 

Fig. 3. Top: Spatial profile of [O iv] emission along the 
major axis of M 82 (symbols) compared to the SWS beam 
profile for a point source (line). Bottom: Radial velocity 
curve for [Oiv] in M 82 (symbols), compared to the CO 
radial velocities as presented by Gotz et al. (1990). A shift 
of -50 km/s has been applied to the CO data. 



sawa et al. |l993j ), but less so for M 82 (Tsuru et al. |1997| ). 
For M 82, definite proof against an AGN origin of [Oiv] 
is provided, however, by the fact that the emitting region 
is spatially resolved and similar in size to the starburst 
region. If it were illuminated by a central source, suffi- 
cient [O iv] could not be produced without exceeding the 
observed relatively low [Ne in]/[Ne n] ratio. 

This constraint is illustrated in Fig. [|, which shows 
predicted line ratios for a simple AGN photoionization 
model with varying ionization parameter computed us- 
ing CLOUDY (Ferland [l996| ). When [Oiv] reaches 1% 
of the low-excitation neon lines, [Nem]/[Neii] is already 
much too high to be consistent with starbursts like M 82 
([Nem]/[Neii]~0.17, Forster-Schreiber et al., in prepara- 
tion). This is a fairly general problem in any photoioni- 
sation scenario, which persists if one adds a small hard 
component to a soft starburst spectrum (as discussed be- 
low in the context of hot stars). The way to circumvent 
this problem - postulate small region(s) with very strong 
[Oiv] but small contribution to the total [Neil] - is not 
viable here since this would imply a central small NLR 
which is inconsistent with the observations of M 82. For 
individual galaxies with [O iv] detections but lacking spa- 
tial information, a weak central AGN remains possible. 
However, we emphasize that the fairly uniform level of 
the [O iv] detections (Fig. ||) requires an unlikely finetun- 
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Fig. 4. Photoionization models for 

[O iv] /( [Neil] +0.44 [Ne in]) (continuous) and 

[Ne in] /[Neil] (dashed) as function of ionization pa- 
rameter in an AGN narrow line region 

n=300, log{U)=-1.5, 39kK Kurucz + 80kK blackbody, stop at T=100K 
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Fig. 5. Photoionization models for 

[O iv] /( [Neil] +0.44 [Ne in]) (continuous) and 

[Ne in] /[Neil] (dashed) in an Hn region excited by 
a 39000 K star, with addition of an 80000 K blackbody 
contributing different fractions of the total Lyman 
continuum luminosity 



ing of AGN and starburst activity to fit our sample as a 
whole. 



4-2. Super-hot stars 

The ionization edge for creation of [O iv] is just beyond 
the He n edge; at higher energies the spectral energy dis- 
tributions of most stars drop precipitously. However, a 
small component of hotter (e.g. Wolf-Rayet) stars might 
provide the necessary high energy photons. We have run 
a photoionization model for an Hn region excited by a 
39000 K main sequence star (represented by a Kurucz 
model atmosphere), plus an additional 80000 K blackbody 
to represent a harder component. The blackbody is an ad- 
hoc choice selected for ease of implementation; however, 
other strong sources of photons beyond 54 eV would give 
similar results. As Fig. ^| shows, this attempt fails to ex- 
plain [O iv] in low-excitation starbursts since the predicted 
[Ne in] /[Neil] ratio (~0.5 to 1) exceeds the observations 
(~ 0.1) when [O iv] reaches 1% of the low-excitation lines. 
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For the high-excitation starbursting dwarfs, such a dis- 
crepancy does not arise, and hot stars remain an option. 

Again, the inconsistency could be alleviated if small 
H II regions with relatively stronger [O iv] emission were 
dispersed in a lower excitation background. In fact, such a 
scenario is qualitatively consistent with the observations, 
as are others with distributed local sources of [O iv] . The 
major reason to consider it unlikely is that we have failed 
up to now to detect [O iv] emission even at a similar level 
in local star forming regions, while we would have to pos- 
tulate regions with stronger emission. The Galactic center, 
which is closest to starburst galaxies in many aspects, still 
shows [Oiv], though even fainter than in the starbursts 
(Lutz et al. |l996b| ). In the massive star forming regions 
W51 IRS2 and 30 Doradus, for which [Ne ill] /[Nell] indi- 
cates high excitation, we were unable to detect [Oiv] at 
a level of <0.01 and <0.005 of [Neii]+0.44[Nein], respec- 
tively (Thornlcy et al., in preparation). 

4-3. Planetary nebulae 

High excitation planetary nebulae are a known source of 
[O iv] emission. A young starburst will, of course, not con- 
tain planetary nebulae and it is easy to show that their in- 
tegrated contribution from the old stellar population is too 
faint. Evolutionary calculations (e.g. Chariot & Bruzual 



1991) show that the contribution of post-AGB stages to 



the bolometric luminosity is less than 1% even in old pop- 
ulations. Making the extreme assumptions that 10% of the 
bolometric luminosity is due to an old population and that 
all PAGB objects are like NGC 7027, one of the highest 
excitation planetary nebulae, we estimate a robust upper 
limit of 10~ 20 Wcm -2 for the [Oiv] emission from plane- 
tary nebulae in M 82, based on [Oiv] flux, luminosity and 
distance of NGC 7027 as given by Shure et al. ( [1983D and 
Beintema et al. ( |1996[ ). 

4-4- Ionizing shocks 

There is ample evidence for ionizing shocks in starburst 
galaxies. Spatially extended, 'Liner'-type optical emission 
lines can be attributed to shocks, and kinematic map- 
ping sometimes provides direct evidence for outflowing 
'superwinds' (Heckman et al. 1990| ). [Oivl column densi- 



ties approaching 10 cm are expected for modest ve- 



locity shocks (100-200 km/s, e.g. Shull & McKee 1979 



Dopita & Sutherland 1996). Assuming postshock values 
of n=1000cm- 3 and T=50000K, we estimate a 25.90/xm 
intensity of ~ 3 x 10 -5 ergs -1 cm -2 sr _1 , equivalent to 
~3x 10~ 20 Wcm~ 2 for the SWS beam. For the assumed 
conditions, the covering factor of such shocks in the star- 
burst region of M 82 would have to be of the order unity. 
At higher shock velocities, the [O iv] column would be in- 
creasingly dominated by material 'at rest' in the photoion- 
ized precursor in the mater ial ahead of the shock front 
(Dopita & Sutherland 1996|) . The shock models predict 



that intensities similar to those estimated for [Oiv] are 
emitted in optical shock tracers like [S ii] 6716/31 A. This 
is fully consistent with optical spectroscopy of M 82 (e.g. 
Gotz ct al. 1990| ). We note that the faint shock emission 
predicted for the [Ne n] and [Ne in] lines will be completely 
dominated by the emission from Hn regions. 

It is instructive to compare the [O iv] results for M 82 
with the SWS observations for RCW 103, a bright super- 
nova remnant interacting with a dense molecular cloud 
(Oliva et al., in preparation). The [Oiv] intensities are 
very similar. The RCW 103 ionic lines are just resolved at 
the SWS spectral resolving power, again similar to M 82. 
Ionizing shocks hence are a plausible origin for the M 82 
[O iv] emission if their total covering factor approaches 
unity in the central starburst region of M 82. 

5. Conclusion 

We have discussed various excitation mechanisms for 
faint [O iv] emission from starburst galaxies. In general, 
starburst-related sources and in particular ionizing shocks 
provide the most plausible explanation. Weak buried 
AGNs may be plausible for individual sources but can be 
ruled out for the best studied case of M 82 whose [O iv] 
emitting region has been spatially resolved. In addition, 
the fairly small scatter in [O iv] versus starburst luminos- 
ity favours a starburst-related origin, since no finetuning 
of two independent mechanisms is required. 
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